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ABSTRACT 


Properties  of  the  hazard  rate  function  are  examined  for 
a  warm  standby  system  composed  of  one  main  component  with 
an  exponentially  distributed  life  length  and  one  standby 
component  with  a  two-stage  exponentially  distributed  life 
length.   Emphasis  is  placed  on  comparing  the  warm  standby 
hazard  rate  function  with  the  hazard  rate  functions  for 
the  parallel  and  cold  standby  systems. 

Graphical  analysis  suggests  certain  hazard  rate  proper- 
ties, some  of  which  are  proved  while  others  are  conjectures. 
Properties  with  possible  applications  examined  include 
initial  and  terminal  values,  monotonicity ,  and  choice  of 
component  for  the  active  role.   Areas  requiring  further 
investigation  are  noted. 
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I.   INTRODUCTION 

Consider  a  system  composed  of  two  independent  compo- 
nents each  of  which  performs  the  same  task.   In  order  for 
the  system  to  complete  its  assigned  mission  either  one 
or  both  components  must  function  throughout  the  mission. 
Each  component  is  assumed  to  have  a  random  life  length. 
A  practical  example  might  be  a  shipboard  steam  generator 
and  an  emergency  diesel  generator. 

The  two-component  system  may  operate  in  one  of  three 
modes : 

1.  Parallel 

2.  Cold  standby 

3.  Warm  standby 

In  the  parallel  mode  both  components  function  at  full  ca- 
pacity simultaneously.   In  the  cold  standby  mode  component 
two  does  not  begin  to  function  until  component  one  fails. 
The  warm  standby  mode  is  the  intermediate  case  where  com- 
ponent two  operates  at  a  fraction  of  its  capacity  until 
unit  one  fails  at  which  time  it  becomes  fully  operational. 

The  hazard  rate  function  is  an  informative  and  intui- 
tively appealing  characterization  of  life  distributions. 
Properties  of  the  cold  standby  hazard  rate  function  were 
studied  by  Barlow,  Marshall,  and  Proschan  [Ref.  1]  when 
the  components  have  increasing  failure  rate  distributions. 
The  hazard  rate  function  for  the  parallel  system  has  been 
studied  by  Esary  and  Proschan  [Ref.  4]  and  by  Bryson 


[Ref.  3]  when  the  components  have  exponential  distributions, 
and  by  Birnbaum,  Esary,  and  Marshall  [Ref.  2]  when  the  com- 
ponents have  increasing  failure  rate  average  distributions. 
The  concern  in  this  thesis  is  to  examine  the  properties 
of  the  hazard  rate  functions  of  the  three  system  modes. 
Component  one  is  assumed  to  have  an  exponential  life  dis- 
tribution.  Component  two  is  assumed  to  have  a  two-stage 
exponential  life  distribution  where  the  first  stage  failure 
rate  is  a  fraction  of  the  second  stage  failure  rate. 


II.   DERIVING  THE  HAZARD  RATE 

A.   THE  SURVIVAL  FUNCTION 

The  survival  function  is  defined  as  the  complement  of 
the  usual  distribution  function.   Let  T>0  represent  the 
random  time  to  failure  of  the  system.   The  survival  func- 
tion is  then  F(t)  =  P(T>t),  t>0. 

For  the  systems  considered  herein,  let  T.  represent  the 
time  to  failure  of  component  i,  for  i  =  1,  2.   Assume  T. 
has  an  exponential  distribution  with  failure  rate  A  .   As- 
sume T   has  the  failure  rate  8A„  for  all  t<T,  and  the  fail- 
ure rate  A_  for  all  t>T1  where  0<6<1.   The  situation  is 
illustrated  in  Figure  1. 

EXP^) 


Component  one 
Component  two 


EXP(0A2)  EXP(A2) 


time 


FIGURE  1.   LIFE  DISTRIBUTION  OF  T 

For  the  system  to  survive  a  mission  of  duration  t, 
either  component  one  must  survive  to  time  t,  or  fail  at 
some  time  s<t  while  component  two  survives  to  time  s  and 
then  on  to  time  t.   Then, 

S/i-\     -Ait  ,  rt   -0Aos   -Ap(t-s)  ,   -Ans  -, 
F(t)  =ex+/e    ze^v      Aie   -1   ds 

o 

-\±t  -A2t  rt.   -(Ai-d-e)  A9)s  ,    ,_^0 

=  e   -1-   +  e   z   /  Aj_e    L  '2'      ds,  t>0. 


The  derivation  proceeds  in  two  cases: 
Case  I.   If  A1=^(l-e)X2,  then 

(A1-(l-e)A2)e-Xlt+Aie"A2t-Aie-(A1+e^)t 
F(t)  =  - = ,   t>0. 

A1-(1-0)A2 
Case  II.  If  X1  =  (1-6) A  ,  then 
F(t)  =  e"Xlt  +  A1te"A2t  ,    t>0. 

i 

This  survival  function  is  applicable  for  all  three  modes 
of  the  system.   That  is,  for  G  =  0,  the  cold  standby  sur- 
vival function  is  Case  I  (A,  ^  A_) . 

,   -Ant    ,   -Ai t 
X-^e      ^   -  A2e   x 

F(t)  = ,    t>0, 

h   -    X2 

Case  II  (Ax  =  A   =  A) . 
F(t)  =  e~At   (1+At) ,    t>0. 
For  6  =1,  the  parallel  mode  survival  function  is 
F(t)  =  e"Xlt  +  e"X2t  -e-(Xl+A2)t  ,    t>0. 


B.   THE  HAZARD  RATE 

The  hazard  rate  function,  r  (t) ,  is  defined  as  the  nega- 
tive derivative  of  the  survival  function  divided  by  the 
survival  function.   That  is,  r(t)  =  f(t)/F(t),  where  f(t)  = 
-dF(t)/dt  is  the  density  function.   For  all  0<G<1  the 
hazard  rate  is 


Case  I. 

Al(A1-(l-e)X2)e"Xlt+A1A2e"X2t-X1(X1+0A2)e"(Xl+eX2)t 


f(t)  = 


r(t)  = 


t>0, 

A1(A1-(l-G)A2)e"Xlt+A1A2e"A2t-A1(A1+6A2)e"(Xl+eX2)t 
(A1-(l-6)A2)e-Xlt+A1e-A2t_Aie-(A1+eA2)t 

t>0. 


Case  II. 


f(t)  =  X±e    Xlt+A1A2te  X2t-A1e  X2t  ,    t>0, 

A1e-Xlt+A1A9te-X2t-A1e-X2t 
r(t)    =  — — ± ,         t>0. 

e-Xlt+A1te-X2t 

For  the  boundary  cases,  the  hazard  rate  function  for 
the  cold  standby  system  (9=0)  is 
Case  I  (A   ^  A  ) . 

-v  i   /^""Aot     "Ait. 
A,  X       (e   *   -  e  x    ) 

r(t)  =  -i-= ,    t>0, 

^  «."^?t   ,   -Alt 
A,e   ^   -A^e  ■L 


Case  II  (A,  =  A   =  A) . 

At 

r(t)  =  ,    t>0. 

1  +  At 

For  the  parallel  system  (9  =  1)  the  hazard  rate  function  is 

A1e-Xlt+A:)e-X2t-(X1+A0)e-(X1+X2)t 

r(t)  =  — - - ,    t>0. 

e-Alt+e-A2t_e-(A1+A2)t 


III.   GRAPHS  OF  THE  HAZARD  RATE 


Norma 


lization  of  the  failure  rates,  A,  and  A  ,  facili- 


tates graphing  of  the  hazard  rate  function.   Let  ex..  = 
A,/(A  +A„)  and  cx„  =  A  /(A,+A  ).   Then  cu  +  a„  =  1.   Also, 
let  t  =  t  (A]L+A  )  .   Then 

Case  I  a,  ^  (1  -  6)a2. 

/      f\       a\  \      "OtlT,        -aoT      .    lQ    >   -  (CXl +0CX?)  T 

a,  (a  -(l-0)a?)e  L   H-a-ja^e   z  -a,  (a,+6a_)e    L        <■' 

i         /n  n\   \  -°hT,    -a?x     -(ai+0ao)T 
(a,- (l-6)a„)e   -*■  +a-,e   z  -a,e    L        z-' 

x>0, 
Case  II  a   =  (1  -  6)a?. 

a,e   x  +anCx~xe   z  -oue   9 
/  \     1        12         1    A       .  _ 
r(x)  =  ,    t>0. 

-air ,     -aoi 
e     +a,ie   ^ 

The  normalization  described  above  is  equivalent  to  the 
change  of  scale  S  =  (A   +A  )T,  for  which 

1             T 
rs(x)  =  : rT( ),    t>0. 

A -I  •   T   A'ry  A  -,    +    A  ^ 

Figures  2.1  through  2.8  are  presented  to  illustrate 
certain  properties  and  conjectures  in  the  next  section. 
The  hazard  rate  and  time  have  been  normalized  in  accord- 
ance with  the  procedure  described  above. 
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FIGURE  2.1.   COLD  STANDBY,  0=0 
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FIGURE  2.2.   WARM  STANDBY,  6  =  .25 
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FIGURE  2.3.   WARM  STANDBY,  6  =  .5 
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FIGURE  2.4.   WARM  STANDBY,  9  =  .75 
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FIGURE    2.5.       PARALLEL,     6=1 
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FIGURE  2.6.   a,  =  .3,  VARYING 
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FIGURE  2.7.   c^  =  .5,  VARYING 


20 


17 


u 


0) 

+J 

U 

rd 

N 

fd 

Q) 
N 
•H 

rH 

rd 

e 
n 
o 

2 


! I — 

10  15 

Normalized   Time,    t 


20 


FIGURE    2.8.       a      =    .7,    VARYING 
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IV.   PROPERTIES  OF  THE  HAZARD  RATE 

The  figures  of  the  previous  section  suggest  that  the 
system  hazard  rate  function  has  certain  properties,  some 
of  which  are  subject  to  easy  proof  and  some  of  which  cur- 
rently stand  as  conjectures. 

A.   THE  INITIAL  HAZARD  RATE 

The  hazard  rate  at  time  t  =  0  is  zero  for  all  three 
modes.   That  is, 
Case  I. 


a,  (a  -  (1-9 ) a0) +a, a  -cu  (a, +9a~) 
r(0)  =  -J— 5 ? 12   11 2_   =  Q/ 

(a  -(l-O)a  )  +  a, -a. 


Case  II. 


r(0)  =  a  +a,a  (0)-a,  =  0. 


B.   TERMINAL  VALUE  OF  THE  HAZARD  RATE 

For  all  three  modes  the  terminal  value  of  the  hazard 
rate  function  is  the  minimum  of  a,  and  a  .   That  is, 

Case  I.   For  a  <a„, 
limr(T)    lim  al(ar(1-G)a2)+ala2e        "al(al+9a2)e 


X  ■*■  oo        T"*"' 


(a1-(l-0)a2)+a1e-(a2-al^-a1e-6a2T 


V 


Case  I.   For  a„<a1, 
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lim  r  (t ) 

X      ->      00 


/         /i    n\       \    -  (on -a?)  7  ,  ,       ,.       .     -  (on -a?  (1-9)  )  7 

,.       a,  (a, - (1-9) a_) e        -      z'    +a,  a„-a,  (a,  +9a~)  e        -1      z  v         'y 
lim      11  2  1211         2 


-f->oo 


(a1-(l-9)a2)e-(al-a2)T+a1-a1e-(al-a2^1-8))T 


=  a2, 


Case  II.   For  all  0<9<1,  a,<a  , 

-(ap-ai)T     -(ao-ai)T 
limr(x)    lim  <»i+<*i"2Te         ~ale    2 

l+a,7e    z   -1' 


C.   MONOTONICITY  OF  THE  HAZARD  RATE 

For  any  value  of  a,  and  a   the  hazard  rate  function  for 
the  cold  standby  system  (9  =  0)  is  monotonically  increasing 
This  is  easily  shown,  and  is  a  special  case  of  a  theorem 
proved  in  [Ref.  1].   If  a,  =  a    ,    the  hazard  rate  of  the 
parallel  system  (9  =  1)  is  monotonically  increasing.   This 
is  also  easily  shown,  and  is  a  special  case  of  a  result 
given  in  [Ref.  4] . 

If  a,  =  a  =  a,  the  hazard  rate  for  the  warm  standby 
system  is  also  monotonically  increasing,  since  then  a,  ^ 
(1  -  9)a2,  that  is  9  ^  1,  and 

r(T)  =  (l+6)g  (1  -  e-9")t    T2Q/ 
1  +  8  -  e-6aT 

from  which 


dr(7)  =  (1+9)  92a2  e  9aT 
dT      (1  +  9  -  e-9aT)2 


/ 


>0,  for  all  7>0. 


This  property  is  illustrated  in  Figure  2.7. 


20 


D.   VARIATION  OF  THE  HAZARD  RATE  WITH  9 


Consider  two  values  of  0,  9,>9„.   One  might  expect  that 


for  fixed  a,  and  a 


rQl(T)>re2  (t)  for  all  x>0. 


Figure  2.7  and  2.8  where  a,  =  .5  and  a,  =  .7  respectively 
suggest  that  this  is  true.   However,  for  a,  =  .3  (Fig.  2.6) 

rei(T)/re2  (x)  for  all  x>0. 

However,  it  appears  that  for  any  fixed  a,  and  a_  and 
any  9>0 


r0(T)>rQ(T)  for  all  t>0. 


This  conjecture  is  suggested  by  Figs.  2.6,  2.7,  and  2 
Moreover,  if  a   =  a   =  a  and  9,>9  ,  then 


rex  (T)>r62  (t)  for  all  x>0, 


since  then 


3r(T)    ae  6aT{  (1  +  9 ) 9ax- (1-e  9aT) } 

,     T  c.  U  , 

89         {1  +  9  -  e~0aT}2 
and  1  -  e"6aT<9ax<  (1+9)  9ax  so  that  8^(t)  >0,  for  all  t>0 

O  V 

E.   INTERCHANGEABLE  COMPONENTS 

Consider  two  components  with  failure  rates  a   and  a 

r(x)a1  =  aA,  c*2  =  aB  =  r(x)  c^-  aR#  a2  =  aA 
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if  0  =  0  or  0  =  1.   For  the  cold  standby  and  parallel  sys- 
tems interchanging  the  two  components  does  not  change  the 
hazard  rate  function  (Figs.  2.1  and  2.5). 

If  a  >a   for  any  O<0<1,  it  appears  that 

r(i)a,  =  a  ,  a   =  a  <r(x)  a,  =  a  ,  a~  =  a   for  all  T>0. 

For  the  warm  standby  mode,  letting  component  one  have  the 
maximum  failure  rate  of  the  two  components  seems  to  yield 
a  lower  system  hazard  rate  for  all  t  than  letting  component 
two  have  the  maximum  failure  rate.   This  conjecture  is 
suggested  by  Figs.  2.2,  2.3,  and  2.4. 

F .   SUMMARY 

The  cold  standby  and  parallel  hazard  rate  functions  have 
been  studied  previously  and  many  of  their  properties  have 
been  proven  and  are  generally  known. 

Because  of  the  more  complex  nature  of  the  warm  standby 
hazard  rate  function  some  of  the  properties  presented  here- 
in are  conjectures  based  on  graphical  analysis.   Further 
investigation  of  these  conjectures  is  needed.   In  addition, 
determining  if  the  warm  standby  hazard  rate  function  has  an 
increasing  hazard  rate  average  for  all  values  of  0,  a. ,  and 
a^   would  be  particularly  interesting.   It  is  known  that 
this  is  the  case  for  the  parallel  system  [Ref.  2]. 
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